Electronic phase separation in iron pnictides 
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A mechanism for electronic phase separation in iron pnictides is proposed. It is based on the com- 
petition between commensurate and incommensurate spin-density-wave phases in a system with an 
imperfect doping-dependent nesting of a multi-sheeted Fermi surface. We model the Fermi surface 
by two elliptical electron pockets and two circular hole pockets. The interaction between a carrier 
in a hole band and a carrier in an electron band leads to the formation of spin-density- wave order. 
The commensurate spin density wave in the parent compound transforms to the incommensurate 
phase when doping is introduced. We show that, for certain parameter values, the uniform state is 
unstable with respect to phase separation. The resulting inhomogeneous state consists of regions of 
commensurate and incommensurate spin- density- wave phases. Our results are in qualitative agree- 
ment with recent observations of incommensurate spin density waves and electronic inhomogeneity 
in iron pnictides. 

PACS numbers: 64.75. Nx, 71.27.+a, 74.70. Xa 



I. INTRODUCTION 

Superconducting iron-based pnictide^ attract consid- 
erable interest not only due to their high critical tempera- 
ture, but also because of the rich physics of their electron 
subsystem. The phase diagram in iron pnictides contains 
areas of superconductivity, spin-density wave (SDW) 
order—, both commensurate^ and incommensurate^, as 
well as the so-called orthorhombic phased. 

In addition, these materials frequently demonstrate 
charge inhomogeneity, exhibiting characteristic features 
of systems with electronic phase separation^"— The ori- 
gin of this phase separation is important for understand- 
ing the mechanisms driving numerous phase transitions 
on the iron pnictides phase diagram. 

The electronic inhomogeneity may be ascribed to sam- 
ple imperfections and thus dismissed. However, another 
explanation to this phenomenon exists. Indeed, it has 
been known that the SDW ground state in a model with 
two spherical Fermi surfaces of unequal radius ^^ can be 
unstable with respect to electronic phase separationi^r— . 
Variations of the latter model have been used to de- 
scribe SDW in chromiun>ii and graphene bilayer with 
AA stackingi^. It is also analogous (but not identical) 
to the commonly used model of iron pnictides where one 
deals, roughly speaking, with elliptical rather than spher- 
ical electron and hole sheets of the Fermi surface. One 
can ask if a similar mechanism could apply to pnictides. 
In this paper, we demonstrate that the answer to this 
question is positive, and thus the charge inhomogene- 
ity can result from a purely electronic mechanism. This 
finding is important for the interpretation of experimen- 
tal data on charge inhomogeneity and for understanding 
the nature of the coexistence of the order parameters in 
iron pnictides. 

For definiteness, we will focus here on iron pnictides. 



although iron chalcogenides also exhibit phase separa- 
tloiiM^ However, the physics related to the Fermi sur- 
face nesting may not be directly applicable to iron chalco- 
genides, especially to those containing alkaline atoms. In- 
deed, some chalcogenides do not have hole pockets, but 
nevertheless exhibit antiferromagnetism with rather high 
Neel temperatures J— ~— Moreover, the electron correla- 
tion effects in chalcogenides seem to be more pronounced 
than in pnictides^. Therefore, the case of iron chalco- 
genides requires a separate consideration. 

This paper is organized as follows. In Sec.[lT]we discuss 
the choice of the model Hamiltonian. Section IIIII deals 
with the study of homogeneous SDW order in the mean- 
field approximation. The instability of the homogeneous 
state is proved in Sec. llVi where the phase diagram is 
constructed as well. The results are discussed in Sec. |Vl 



II. MODEL 

A. Kinetic energy 

Unlike cuprates, which are believed to be in the strong 
electron-electron interaction regime, iron pnictides may 
be described by a weak-interaction model (see, e.g., dis- 
cussion in Sec. HI A of Ref. 12 ll). In this approach the 
shape of the Fermi surface and the value of the Fermi 
velocity are the only relevant single-electron band pa- 
rameters. The commonly used form of the Fermi surface 
for iron-based superconductors is shown schematically in 
Fig. [TJa)^— Here, we plot our Fermi surface within 
the first Brillouin zone, which corresponds to the square 
lattice of iron atoms, with one Fe atom per unit cell 
and the lattice constant a. Two quasi- two-dimensional 
nearly-circular hole pockets centered at the r(0, 0) point, 
and two elliptical-shaped electron pockets centered at 
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Figure 1: (Color online) (a) The Fermi surface of iron pnic- 
tides. Two circular hole pockets are centered around the 
r(0, 0) point, and two elliptic electron pockets are located 
near the (0,7r/a) and (7r/a,0) points. According to Ref. 22, 
only one hole and one electron band are involved in the for- 
mation of the SDW with a nesting vector Q = Qo + q, where 
Qo = (^r/a, 0), while two others remain untouched, giving rise 
to the metallic state for both the parent and the doped com- 
pound, (b) Fermi surface nesting at different doping levels 
for a > 0: the electron pocket is shown shifted by the vector 
— Q. For the parent compound q = 0, while for the electron 
(hole) doping x > (x < 0), the vector q is parallel to the x 
axis (y axis), and qx{x) = qy{—x) (see text). 



the M(0,7r/a) and M(7r/a, 0) points are taken into ac- 
count. Less important parts of the Fermi surface are 
ignored. We will study a two-dimensional model, assum- 
ing that the interlayer electron tunneling and interactions 
are smaller than those inside the two-dimensional layers. 
Thus, the Hamiltonian of the model has the form 

H = Ho + Hinu (1) 

where the kinetic energy term Hq is given by 

H0 = Y1 4k«L.«k.. + E ^tAsaKs. ■ (2) 

kscr kscr 

In this equation, a^^^, a^^^ (^kscr' ^kscr) ^^^ ^^^ creation 



and annihilation operators for electrons in the hole-like 
(electron-like) bands s = 1,2 with the spectra e^^ (^sk)- 
To simplify our formalism, we assume that near the 
Fermi level the bands have quadratic dispersions. For 
the central circular hole-like bands we have {h = 1) 
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where kp is the Fermi momentum, /i is the chemical po- 
tential, v^ and Vp are the Fermi velocities for the hole 
bands. The energy shift Ae in Eq. ([3|) defines the dif- 
ference in radii of the two hole pockets: for Ae > the 
radius of the second hole pocket is smaller than that of 
the first one. 

As mentioned above, the electron components of the 
Fermi surface are elliptic. For these, the dispersion near 
the Fermi surface is given by the following relations: 



^Ik+Qo — ^K^xiky)^ ^2k+Q^ — ^{ky^kx)^ 



(5) 



where the function e is equal to 



^\f^xi f^y) — 



'J-pi^F 
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^^{^l-kl)-^. (6) 
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The centers of the elliptic bands are Qo = (7r/a, 0) and 
Qo = (0,7r/a), and v^ is the Fermi velocity for the elec- 
tron bands averaged over the Fermi surface. Note that 
in general v^ ^ Vp ^ Vp. The chemical potential fi of 
the undoped compound is zero. The parameter a defines 
the ellipticity of the electron pockets. For the Fermi sur- 
face structure corresponding to Fig. [IJa), a is positive. 
In this case, the major axes of the ellipses are directed 
toward the F point. For a < 0, the ellipses are rotated 
by 90° around their centers. 



B. Interaction Hamiltonian 

Due to the multi-sheeted structure of the Fermi sur- 
face, the interaction Hamiltonian i^int, in general, must 
include a number of terms describing interactions be- 
tween carriers in different bands. However, since we are 
interested in the SDW order, most of these terms may 
be omitted for they do not contribute to the SDW phase 
transitionii. For this reason we ignore the electron- 
electron and hole- hole interactions. 

It is known^ that magnetic order in pnictides is 
"stripy". That is, the value of local magnetic moment 
oscillates along one of the crystal axes, remaining con- 
stant along the other axis. Such magnetic structure cor- 
responds to the order parameter, which couples one of 
the hole bands (for definiteness, assume that it is the 
hole band e^) with one of the electron bands (assume 
that it is ef)^^. The order parameter, which couples e^ 
with both £\ and 82-, is unstable. ^^ Consequently, only 



the interaction between s^ and el is crucial for the stabi- 
Hzation of the ordered phase. All other interaction terms 
will be discarded. We will re-examine this simplification 
in subsection IIIIBI [see discussion after Eq. (|36]) ]. 

Keeping these considerations in mind, we split the 
model Hamiltonian into two parts, magnetic, i^m, and 
nonmagnetic ("reservoir"), H^: 



H = 



Hm + Hr , 



(7) 
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Here M is the number of Fe atoms in a layer, Vi > is the 
coupling constant characterizing the Coulomb interaction 
between the bands e^^ and £\^. Below we will refer to 
these bands as magnetic bands. The bands £2k ^^^ ^2k 
will be called nonmagnetic, since these do not contribute 
to the magnetic order parameter. 

Model d?]) is a generalization of the Rice model, pro- 
posed in Ref. [ll| for the description of the incommen- 
surate SDW order in chromium. The Hamiltonian of 
Ref. [ll| has two bands, which participate in the mag- 
netic transition, and a "reservoir" (nonmagnetic bands 
corresponding to i^r, where the subscript 'r' stand for 
"reservoir"). Unlike Eq. ([7j), the bands responsible for 
magnetic ordering in the Rice model have a spherical 
Fermi surface; therefore, at certain filling, the nesting is 
perfect. 



III. INCOMMENSURATE SDW ORDER 

A. Mean-field equations 

We now consider the Hamiltonian ([7j) in the mean-field 
approximation. As we pointed out above, the Coulomb 
interaction in iron pnictides is weak, therefore, we assume 
below that 

Vi/SF^l, (10) 

where the Fermi energy sp is defined as 



sf 



Vr + V'r 



■kp = vpkp ' 



(11) 



The weak-coupling condition Eq. (p!Q]) guarantees the ap- 
plicability of the mean-field approximation. 

We will study the stability of the following SDW order 
parameter: 



Vi^ 
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(12) 



where the nesting vector Q is equal to 
Q = Qo + q. 



(13) 



When q = 0, our SDW is commensurate, whereas if q is 
small, but non-zero, it is incommensurate. Other types 
of order parameter will be discussed in Sec. [Vl 

The magnetization corresponding to the SDW order 
parameter A lies in xy plane. For the commensurate 
SDW case, in real space the "stripy" order is observed. 
Namely, the magnetization direction remains constant 
when one moves along the direction normal to Qo- How- 
ever, when one moves parallel to Qo, the magnetization 
reverses its direction from one iron atom to the next iron 
atom. For incommensurate SDW, this "stripy" pattern 
slowly rotates in the xy plane: the local rotation angle 
^(R) at the point R is equal to (qR) [see panel (c) and 
(d) in Fig. [2]. 

In the mean-field approximation, the magnetic Hamil- 
tonian Eq. (jHI takes the form 



H^ 



MF 



5Z [^lk-q4aCk,+ ^?k+Qo4a^ka 
kcr 

- A (cLdk-. + 4.Ck-.) + AVV^i] , (14) 



where we introduce the new operators cj^^ = CLl-_qia^ 
dl.^ = ^k^Q 1^, and A is assumed to be real. The Hamil- 
tonian Eq. (fT^ can be easily diagonalized. The quasipar- 
ticle energies are 



E, 



(1,2) 
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The grand potential Qm (per one Fe atom) corresponding 
to H^^ is {kB = 1) 
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where vq is the volume of the unit cell, and the integra- 
tion is performed over the 2D Brillouin zone. The SDW 
gap A and the nesting vector Q are found from the min- 
imization of r^m 



dOrn 

dA 



0, 
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(17) 



It will be shown below that in the weak-coupling limit 
Eq. (fTQ|) . the SDW order can exist only if the deviation 
from the perfect nesting is small, that is: 



lal <C1. 



(18) 



In addition, the gap A is small compared to Sf^ and the 
deviation q of the nesting vector Q from the commensu- 
rate value Qo is also small 



IQ-Qol ^ A/vF<^kF. 



(19) 



Restricting ourselves to the limit of zero temperature, 
we can write the first of Eqs. ([T7j) in the following form 






(20) 



A2 + -(£?k+Q„-^?k-q) 



where 9(x) is the step function. 

When the lower band E^. is filled, while the upper 

(2) 

band E^ is empty, the gap attains its maximum value 

Aq. In this case, e{E^^) = Q{-E^^) = for any k, and 
Eq. ([2Q|) becomes 



1 = ^ IdE- 



P{E) 



2 J .MTW' 

where the generalized density of states is defined as 



(21) 



where, 



^^o(p,V^) = Px cos (f^py simp - -cos2(f, (27) 



a 



a ksf 



(28) 



Ao ' 
and we introduce the following dimensionless quantities 



2^w^4 

Vf + ^F 



^ = -^' i^ = ^^, P = ^^ ♦ (29) 
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Transforming similarly the second of Eqs. (p!7|) . we 
obtain the equation for the nesting vector Q = Qo + 

2Aop/vf^, 
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(22) 



Evaluating p{E) we set q = 0, since taking into account 
values of the order of |q| ~ A/vf gives only second-order 
corrections in Eq. ([2T]) . 

If \E\ <C ^F, the function p{E) can be calculated ex- 
plicitly using Eqs. (J3]) and ([5]) for the band spectra ^ik+QQ 
and ej^j^. As a result, near the Fermi surface we obtain 



piE) « m 
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\E\ <£F. 



(23) 



When the energy E is of the order of the band width, 
p{E) vanishes. This makes the integral in Eq. ([2T]) con- 
vergent, and one can derive the usual BCS-like expression 
for the gap: 



Ao ~ Sf exp 



27r£F 



(1) 



(24) 
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If the sample is doped, then E^^ > or E^ < 
for some range of k, and the equation for the band gap 
becomes 
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Substituting Eqs. ([3|) and ([5|) into Eq. ([25]) . and taking 
into account Eqs. (fTSj) and ([19]), after straightforward 
algebra we derive the equation for the band gap in the 
form 



27r 




z^o(p,V^) 



(26) 



Re\/{uo{p,(p) -u) -(^2. 



(30) 



Equations ([26]) and ([30]) determine the SDW band gap 
A and the nesting vector Q as functions of p. However, 
experiments are performed at fixed doping, not chemical 
potential. Thus, we have to relate the electron density 
and /i. The total number of electrons per iron atom, n{p), 
is the sum of the number of electrons in the nonmagnetic 
and the magnetic bands: n{jii) = nr(/i) + ^m(M), where 



2r(/i) 



rimip) 



M 



E[©(-^2k) + e(-e^i,)] 

k 
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The doping is defined by 

x{p) = n{p) — n(0). 



(31) 



(32) 



Performing calculations similar to those described above, 
we obtain for the doping level 



Xo 



27r 

/^ sgn(i/o(p,9^) -1^) X 



^cJ{iyo{p,ip)-u) -S^, 



where 



Xo 
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TTK Sf ' 






2xlv^pv'p 



(33) 



(34) 



The first (second) term in Eq. ([33]) is the nonmagnetic 
(magnetic) contribution x^ {x^) to the total doping x. 
Equations ([26]) . ([30]) . and ([33]) form a closed system of 
equations for the self-consistent determination of A(x), 
Q(x), and p{x). 



B. Results: Homogeneous state 



Our numerical analysis reveals that Eq. ([26]) has no 
solutions if |a| > 2.0. Thus, the equality |a| = 2.0 deter- 
mines the critical value of the ellipticity parameter olc'. 
for a given coupling constant V\^ the SDW ordering oc- 
curs only if 



\a\ < an 



2Ao 

KSf 



(35) 



Conversely, this condition Eq. (|35]) may be re-formulated 
as a requirement on the interaction strength: for given 
band parameters (a, x, etc.) the SDW is stable only if 



Vi>Vc = 



27r€F 



vokjp In 



(36) 
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Equation ()36|) allows us to re-examine the assumption 
we made in subsection III B[ When constructing the inter- 
action Hamiltonian, we discarded several terms including 
the interaction between electrons in ef 2 ^^^ holes in the 
second hole band 62 • When mapping the phase diagram, 
Eq. (|36|) is very useful since it assures that a sufficiently 
weak interaction does not affect the phase diagram: to 
induce a new ordered phase, the corresponding coupling 
constant must exceed Vc- Since the ordering of the sec- 
ond hole band is not observed experimentally, we con- 
clude that the pertinent interaction is too weak, and, 
consequently, may be neglected. 

It is seen from Eqs. ([29]) and ([34]) that the parameter 
K only renormalizes the dimensionless quantities u^ p, 
and xq. Thus, at fixed a and r, it is enough to know the 
functions A(x), q(x), and fi{x) for x = 1. For x 7^ 1 
these functions are found by simple rescaling as A(x>^), 
q(xx)/x, and KjiixK). For this reason, below we study 
only the x = 1 case. Note that, typically, for a real 
material all Fermi velocities are of the same order, thus 
H ^ 1. 

The computed functions A(x) and q{x) are shown in 
Fig. [2]^a) and (b) for r = x = 1 and two different values 
of a. At low doping, \x\ < x^, where the critical doping 
xl depends on model parameters, all the extra charge 
goes to the nonmagnetic bands. As a result, the order 
parameter is independent of x, A(x) = Aq, q(x) = 0, 
and the chemical potential increases linearly with x. 



li{x) 



AqX 

rxo ' 



(37) 



In such a regime, the system shows commensurate SDW 
order. 



When 



the chemical potential touches the 



bottom (top) of the upper (lower) magnetic band. Upon 
further doping, electrons (holes) appear in the band E^ ^ 

(E^. ^). The SDW order becomes incommensurate^ q 7^ 
0. The order of the transition into the incommensurate 
state depends on the parameters a and r. If r < ri = 
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Figure 2: (Color online) Dependence of the normalized order 
parameter A/Ao and wave vector |q| on the doping x, for r = 
K — 1. (a) Here a — 0.1. The transition from commensurate 
to incommensurate SDW state is of first order: both A and 
|q| abruptly change at x = xj. The dashed curves show the 
functions A(x) and |q(x)| corresponding to the metastable 
state, (b) Here a = 0.6. The transition from commensurate 
to incommensurate SDW state is now of second order. Panels 
(c) and (d) show schematics of the incommensurate SDW 
spin structure. If a > 0, then panel (c) corresponds to hole 
doping (x < 0) and (d) to electron doping (x > 0). If a < 0, 
then panel (c) corresponds to electron doping, while panel (d) 
corresponds to hole doping. 



0.38, the transition is of second order for any a, while for 
larger r it becomes of first order if \a\ < ai{r). In the 
case of first-order transitions, both A and |q| abruptly 
change at x = x^, see Fig. [2](a). 

Electron doping {x > 0) is in many respects similar to 
hole doping {x < 0). Indeed, the functions A(x), /i(x), 
and |q(x)| obey particle- hole symmetry relations: 



A{-x) = A(x), 
/i(-x) = -/i(x), 

|q(-x)| = |q(x)|. 



(38) 
(39) 
(40) 



In addition, these functions are independent of the sign of 
a. However, the direction of the vector q(x) depends on 



the sign of a and on the type of doping. When a > 0, the 
vector q(x) is parahel to the x axis for electron doping 
and to the y axis for hole doping. The direction of q(x) is 
reversed if a < 0. Thus, the nesting vector Q = {Qx-, Qy) 
is given by the following formulas 



a kvf 

^ _ 2Aop(|x|) 

1 ^y — 

a kvf 



= 0, if ax > 0, (41) 
, if ax < 0. (42) 



How does the incommensurate SDW described by these 
equations look like in real space? Let us introduce the 
integer- valued vector n = (n^m). Then the position of 
a given iron atom is equal to Rn = an. For the order 
parameter A, Eq. (p!2]) . the SDW magnetization vector 
lies in the xy plane: S^ = 0. The in-plane components 
can then be expressed as: 

Sn = 5*0 ( cos(aQ • n), sin(aQ • n) j , (43) 

where So = A/Fi. Thus, we obtain 

g ^Q f (-l)''(cos0n, sin^n) , ax > , .^^. 
" ^1 (— l)^(cos(/)7Ti , sin^TTi) , ax < , ^ ^ 

where (j) = a|q|. These spin configurations are schemat- 
ically shown in Fig. [2] (c) and (d) for hole and electron 
doping, respectively. For positive (negative) a, panel (c) 
corresponds to x > (x < 0), and panel (b) corresponds 
to X < (x > 0). 

IV. PHASE SEPARATION 

In the previous section we assumed that the ground 
state of the model is homogeneous. Here we demonstrate 
that there is a part of the phase diagram where homoge- 
neous states are not stable, and the true ground state is 
phase separated. To detect this instability, the chemical 
potential must be calculated. 

The computed dependence of the chemical potential /i 
on doping x is shown in Fig. [3] for three different values 
of a and r = x = 1 for electron doping {x > 0). The 
dependence /i(x) for hole doping can be easily determined 
using particle- hole symmetry, as discussed in the previous 
section. 

The function /i(x) shown on Fig. [3] demonstrates sev- 
eral peculiar features in the vicinity of the commensurate- 
incommensurate phase transition point x^. Specifically, 
for low a, the function /i(x) is non-monotonic and multi- 
valued near x ^ xl [see Fig. [31(a)]. At higher a, the mul- 
tivaluedness disappears; however, the non-monotonicity 
remains [see Fig. [31(b)]. This vanishes at even higher val- 
ues of a [as in Fig. [31(c)]. 

We now observe that in Fig. [S^a) and Fig. [S^b) there 
are finite ranges of doping, where the chemical potential 
decreases as the doping increases. This means that the 



compressibility of the electronic system is negative and 
the homogeneous state is unstable with respect to sepa- 
ration into two phase8^. The phase transitions between 
homogeneous commensurate and incommensurate SDW 
phases, which we described in subsection IIII B\ will be 
masked, at least partially, by the phase separation. 

In the separated state there exist two phases, phase 1 
and phase 2, with electron density xi{< x^) and X2{> 
xl)^ and with the volume fractions pi and p2 satisfying 
the conditions pi + P2 = 1 and xipi + X2P2 = x. As one 
can see from Fig. [3l (a,b), the phase 1 is the commensu- 
rate and the phase 2 is the incommensurate SDW state. 
The concentrations xi and X2 are found from the equa- 
tions^-^ fi{xi) = fi{x2) = /io and Qi = 1^2, where 1^1,2 
are grand potentials in the phases 1 and 2. The latter 
condition means the equality of two shaded areas shown 
in Fig. [3l(a,b) (the so-called Maxwell construction). 

The range of doping x, where the phase separation 
exists, is largest when a = 0. In this case our model 
is identical to the two-dimensional Rice model^^, for 
which the presence of the phase separation was shown 
in Refs. 12.13. The range of the phase separation, 
Xi < X < X2, shrinks if |q;| increases, and disappears at 
the critical value ac — 1.15. The phase separated state 
does not exist for |a| > ac- 

The results obtained are summarized in the phase di- 
agram in the (x, a) plane shown in Fig. [H This phase 
diagram is calculated for r = x = 1. It remains qual- 
itatively the same if r 7^ 0. If the nonmagnetic bands 
are absent, r = 0, the homogeneous commensurate SDW 
phase exists only when x = 0. Consequently, the elec- 
tronic concentration xi in the phase 1 is zero in the phase 
separated state for any a. 

Finally, we would like to draw the attention of the 
reader to an unexpected feature of the phase diagram in 
Fig. [H Namely, the second order incommensurate SDW- 
to-paramagnet transition line bends to the right, while 
all other transition lines on the phase diagram bend to 
the left. In other words, the value of doping where the 
transition from paramagnetic to incommensurate SDW 
phase occurs increases when the ellipticity parameter a 
increases. Interpreting this feature, however, one must 
keep in mind that the nesting quality is controlled not 
only by a but also by /i (or x). The parameter a con- 
trols the shape of the electron pocket, while fi (or doping) 
controls the relative areas of electron and hole pockets. 
Therefore, at larger doping levels, when the areas of the 
electron and hole pockets differ substantially, the elec- 
tron and hole pockets could be better nested for larger 
ellipticity of the electron pocket. 



V. DISCUSSION 

The most important result obtained here is the pre- 
diction of electronic phase separation in some range of 
doping. The separated phase consists of a mixture of 
commensurate and incommensurate SDW phases with 
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Figure 3: (Color online) Chemical potential \i(x) calculated 
at q; = 0.1 (a), a — 0.5 (b), and a — \fl (c). Here, r — k —\. 
(a) and (b): the homogeneous state is unstable toward phase 
separation if xi < x < x^. The dashed (red) curve corre- 
sponds to the /io found using the Maxwell construction. The 
shaded areas above and below /xo are equal to each other, (c): 
\i{x) monotonically increases with x, no phase separation ap- 
pears. The homogeneous commensurate (q = 0) and incom- 
mensurate (q / 0) SDW, paramagnetic (PM), and inhomoge- 
neous commensurate-incommensurate SDW (PS) states are 
separated by vertical dotted lines. 



different electronic concentrations. The phase separa- 
tion in iron-based superconductors was observed in sev- 
eral experiments^"—. For example, the inhomogeneous 
state with a commensurate antiferromagnetic and non- 
magnetic domains with characteristic sizes ~ 65 nm was 
observed in the hole doped Bai_a;Ka;Fe2As2 compound.^ 
Our theory predicts that the second phase is the incom- 
mensurate SDW rather than a nonmagnetic one. How- 
ever, the proposed mechanism of phase separation can, 
in general, be consistent with the observations reported 
in Ref. |6|. We found that the thermodynamic potentials 
of the incommensurate SDW and the metastable param- 
agnetic phases are very close to each other in the doping 
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Figure 4: (Color online) The phase diagram of the model d?)) 
in the (x,a) plane; for r — >c — \ and a, x > 0. It is 
symmetric with respect to the replacement x ^ —x and/or 
a ^ —a. The boundary between incommensurate SDW and 
paramagnetic (PM) states shown by the dashed (red) curve 
corresponds to the first order and by the solid (red) curve to 
the second order phase transition. Solid (blue) lines indicate 
the boundaries of the phase separated (PS) state. The solid 
(black) curve (second order transition) and dashed (black) 
curve (first order transition) show the boundaries between 
commensurate (q = 0) and incommensurate (q 7^ 0) homo- 
geneous SDW phases. Note that phase separation partially 
masks the transition line between the homogeneous SDW 
states. 



range x > X2' The incommensurate SDW phase can be 
destroyed by an additional reason not taken into account 
in our model, e.g., by disorder. In this case the phase 
separation might occur between the commensurate SDW 
and the paramagnetic phases. In addition, the incom- 
mensurate SDW order parameter may be difficult to de- 
tect due to its weakness. 

The geometry of the emergent inhomogeneities and 
their characteristic sizes are beyond the scope of the 
present study. The electron concentrations in sepa- 
rated phases are different; hence, the inhomogeneities 
are charged and one should take into consideration the 
electrostatic contribution to the total energy. The char- 
acteristic sizes of the inhomogeneities are controlled by 
the interplay between the long-range Coulomb interac- 
tion and the energy of the surface between the phases. 
In the simplest case, the structure of the inhomogeneous 
state corresponds to the droplets of one phase embedded 
in the matrix of another phase. However, depending on 
the properties of the system, other geometries are possi- 
ble, such as alternating layers of different phases, stripes, 
etc^— When disorder is present in the sample, it also 
affects the structure of the inhomogeneous phase. 

In our study we chose the order parameter A, Eq. (fT2]) . 
In the literature, a different order parameter is also 
discussed J^i^ This order parameter is not homogeneous: 
SDW gap experiences periodic modulations in real space. 



When the state with such an order parameter is doped, 
the extra charges go to the places where the gap locally 
vanishes. The latter type of order may be more stable 
than Eq. ([12]). However, for several reasons we decided 
to avoid such an option. Specifically, analytical calcula- 
tions with inhomogeneous gap become very complicated. 
Further, the choice of the order parameter affects some- 
what the phase diagram and, in particular, the region 
of the phase separation; however, the inhomogeneous re- 
gion does not disappear. In addition, we must remember 
that the relative stability of different order parameters is 
likely a non-universal quantity, which depends on a vari- 
ety of microscopic parameters (e.g., details of the band 
structure, interaction, disorder). Therefore, the type of 
suitable order parameter cannot be deduced without in- 
put from experiments. 

It follows from the phase diagram in Fig. [4] that the 
smaller effective ellipticity parameter a, the larger the 
region of the phase separation. The SDW gap Aq in- 
creases with increasing coupling potential Vi. Therefore, 
a stronger coupling is more favorable for phase separation 
since a ex a/Ao [see Eq. (|28|) ]. 

Our mean-field approach is not applicable in the limit 
of strong coupling, where the use of Hubbard-like models 
is more appropriate. The phase separation in such mod- 
els is a common phenomenon^—. Thus, we expect that 
the existence of phase separation in our model is not an 
exclusive feature of the weak-coupling regime. 

The model studied here predicts an incommensurate 
SDW phase for both electron and hole doping in the 
range of dopant concentrations scaled by the parameter 
xo, Eq. ([34|) . From the latter equation we conclude that 
xo <; 1 if the weak-coupling condition, Eq. (fTQ|) . is met. 
Thus, we predict the existence of the incommensurate 
SDW order for small doping, x <C 1. 

The incommensurate SDW phase is characterized by 
the vector q = Q — Qo- Our calculations show that there 
are only two possible equilibrium directions of q: it can 
be either parallel, or perpendicular to Qo, depending on 
the type of doping and the sign of the ellipticity param- 
eter a [see Eqs. (|4T]) and (|42]) ]. It is clear, however, that 
both the magnitude and direction of q is sensitive to the 
shape of the Fermi surface, which is simplified in the 
present calculations. In real materials, the shape of the 
hole pockets deviates from perfect circles. Moreover, the 
spectrum of the hole and electron bands depend on the 
transverse momentum, /c^, which is completely neglected 
in the model considered here. In particular, the electron 
bands have corrugated structure^i^, that is, the major 
axis of the ellipses are rotated by 90° when kz varies from 



to TT. Nevertheless, we believe that the model studied 
here captures the main features of the incommensurate 
SDW state in pnictides. 

Indeed, the observation of the incommensurate SDW 
phase with q perpendicular to Qo in the electron doped 
Ba(Fei_a^Coa;)2As2 was recently reported in Ref. 0. The 
measured value of |q| was about 0.02-0.03 for the con- 
centration range 0.056 < x < 0.06. Such values of x and 
q, as well as the direction of q (corresponding to a < 0), 
are consistent with our theory. 

The short-range incommensurate SDW phase was de- 
scribed in Ref . l48l for the electron doped BaFe2-a^NiajAs2. 
In general, one must be cautious interpreting this exper- 
iment using our model, since the observed short-range 
SDW correlations do not correspond to our long-range 
SDW order. However, we note that the measured vector 
q was perpendicular to Qo, which is also consistent with 
a negative a. 

We found no experimental work reporting the observa- 
tion of the incommensurate SDW phase in the hole-doped 
pnictides. However, double-peaked spin fluctuations at 
Q = Qo ± q have been observed^S*^ in the hole-doped 
Bai_a;Ka;Fe2As2 compound in a wide doping range. In 
these measurements^!^ the vector q was found to be 
parallel to Qo, which, again, corresponds to a < in our 
model. 

In conclusion, we study a model with imperfect nesting 
of the Fermi surface suitable for the description of iron 
pnictides. We show that the incommensurate SDW phase 
arises in the system at a finite doping level. We demon- 
strate that two spin configurations can arise in the sys- 
tem depending on the model parameters. It was shown 
that the homogeneous state is unstable toward phase sep- 
aration into commensurate and incommensurate SDW 
phases in a specific doping range. These results are in 
qualitative agreement with the recent experimental ob- 
servation of the incommensurate SDW order and phase 
separation in doped superconducting pnictides. 
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